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We derive a nonlinear Schrodinger equation with a radical term, ±-\/l — \V\'^ , as an asymptotic 
model of the resonantly absorbing Bragg reflector (RABR), i.e., a periodic set of thin layers of two- 
level atoms, resonantly interacting with the electromagnetic field and inducing the Bragg reflection. 
A family of bright solitons is found, which splits into stable and unstable parts, exactly obeying 
the Vakhitov-Kolokolov criterion. The soliton with the largest amplitude, (IV^Djn^x = 1, is a "quasi- 
peakon", i.e., a solution with a discontinuity of the third derivative at the center. Families of 
exact cnoidal waves, built as periodic chains of quasi-peakons, are found too. The ultimate solution 
belonging to the family of dark solitons, with the background level V = 1, is a dark compacton. 
Those bright solitons which are unstable destroy themselves (if perturbed) attaining the critical 
amplitude, \V\ =1. The dynamics of the wave field around this critical point is studied analytically, 
revealing a switch of the system into an unstable phase, in terms of the RABR model. Collisions 
between bright solitons are investigated too. The collisions between fast solitons are quasi-elastic, 
while slowly moving ones merge into breathers, which may persist or perish (in the latter case, also 
by attaining \V\ — 1). 

PACS numbers: PACS numbers: 42.25.Bs, 42.50.Gy, 78.66.-w, 42.65.Tg. 



I. INTRODUCTION AND THE MODEL 



The interplay between the resonant reflection of light on Bragg gratings (BGs) and resonant interaction of light with 
nanolayers of two-level atoms (of width < 100 nm, which is much smaller than the wavelength of hght), or with similar 
active elements, deposited at reflecting layers of which the BG is built, gives rise to artificial optical media in the 
form of resonantly absorbing Bragg reflectors (RABRs). They are promising for fundamental studies and applications 
including the storage of slow-light pulses [1], negative reflection 0], and periodically amphfying BGs Q. The 
currently available nanofabrication techniques make the creation of RABRs with required properties quite feasible . 
In particular, the interplay of the resonant nonlinearity (which gives rise to the self-induced transparency in uniform 
media [31) with the bandgap spectrum induced by the BG rnay give rise to peculiar species of temporal solitons in 
RABRs (see Ref. [l^l for a review, and more recent works [ll| - [l3j ). These solitons belong to the class of gap solitons, 
whose propagation constant falls into the underlying bandgap. Gap solitons in fiber Bragg gratings with the uniform 
Kerr nonlinearity have been a subject of intensive theoretical studies and have been created experimentally, in 
the form of moving Bragg solitons 

Using the rotating-wavc approximation, the Maxwell-Bloch equations governing the transmission of light in the 
RABR can be reduced to the system of equations for the scaled slowly varying variables, viz., the amplitude of the 
electromagnetic field, S^, polarization of the medium, P, and population inversion w of two- level atoms [3l.[l0j: 

(S+),, - (S+)^^ = + 2lTjP + 2Pr, (1) 

Pr = -iSP + S+w, Wr = -Re i^+P*) , (2) 

where r and ^ are the scaled time and coordinate (another amplitude, E_, is governed by a detached equation). 
Coefficient S, which may be positive or negative, measures the detuning of the transition frequency of the two-level 
atoms from the carrier frequency of the electromagnetic field, and 77, which is defined to be positive, is the scaled BG 
reflectivity. 

A straightforward corollary of Eqs. ([2]) is 9 ^1^*1^ + /9r = 0, i.e., \P\^ + = const. The normalization may 
be set by fixing const = 1, hence w can be eliminated in favor of P 0]: 

w = ±v/l-|PP. (3) 

The stable situation is determined by the condition that, in the absence of the polarization (P = 0), the atomic 
population must be uninverted {w = —1), hence the lower sign must be chosen in Eq. ([3]). The opposite situation is 
possible too, but it is unstable, corresponding to an inverted population in the absence of the field. 

Thus, assuming that the stable branch of square root ([3]), is used, w = — i/l — |Pp, one arrives at the system of 
two equations Q, 



(S+)^^ - (S+),, - -r;2s+ + 2z(r; - 6)P - 2^1 - |PP S+, (4) 
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Pr = -i5P- ^Jl-\PV S+ . (5) 

Looking for bright-soliton solutions to Eqs. ^ and ([S]) as S+ = exp (— iwr) 5(C), f = i exp (— iwr) 7^(0, solutions 
for real functions 5(C) and P(C) were found in an implicit analytical form in Ref. In the general case, these 
solitons fall into two distinct bandgaps produced by the linearized version of the system. Dark-soliton solutions were 
also studied in Ref. Q ■ 

The present paper is focused on subfamilies of solitons residing near the edge of one of the gaps, oj = 5. First, we aim 
to shown that, in an asymptotic approximation valid in this case, Eqs. (j4]) and (O reduce to a single nonpolynomial 
Schrodinger equation (NPSE), with the nonlinear term in the form of a radical. To this end, solutions are looked for 
as 



e 



^Q-iy(C, r), P = -le-^'^ [y(C, t) + i?(C, r)] , (6) 



where Q = {q + 5) /2, y (C,t) is assumed to be a slowly varying function of time, in comparison with exp(— zi5t), 
and i? is a small correction to V required by the self-consistent derivation. Actually, the slow time dependence in V 
accounts for a small deviation of the full frequency from uj ~ 5] in other words, for the purpose of the asymptotic 
analysis 5 and 77 may be considered as large parameters, which corresponds to a far-detuned strongly reflecting RABR. 
Then, it is a simple exercise to demonstrate that the self-consistent asymptotic approximation corresponds to 



R = &-^[iVr-Q-^^/l^\VVV) (7) 
in Eq. ([S]), and the asymptotic NPSE takes the following form: 



iVt + Va - e^JT^\V?y - 0, (8) 



where the rescaled time \s t = 5 {rf + 5^) ^ r, and e = 2r]/5. An additional obvious rescaling of t and C allows one 
to fix e = ±1 for (5 ^ 0, which is adopted below. It is easy to see that Eq. ([5]) gives rise to bright- and dark-soliton 
solutions for e = +1 and —1, respectively, i.e., for positive and negative values of the mismatch, (5 > and 5 < Q. 
Equation ([8|) conserves three dynamical invariants, viz., the momentum, Hamiltonian, and norm, 

/+00 
mOl'dC- (9) 
-00 

Stationary solutions to Eq. ([8]) (in particular, solitons) arc looked for in the usual form, 

F(C,i)=e"*'^*W^(C), (10) 

where real function W{Q satisfies equation 

^ = -,«' + ../T^H.-.-5|l. (11) 

with effective potential t/cff = (l/2)xW^^ + (£/3) (l ~ W'^)^^^ . Substituting such stationary solutions back into Eqs. ([6]) 
and ([7]), one can reproduce the respective solutions in the framework of the underlying RABR model. In particular, it 
is worthy to note that real stationary solutions W{C) correspond, according to Eqs. ^ and ([7]), to complex stationary 
states of the material polarization, i.e., those with an intrinsic chirp. 

Another noteworthy finding reported below, which directly pertains to the relation between solutions to the NPSE 
([HI) and solutions of the underlying system of Eqs. ^ and ([2]), is that when, in the course of the dynamical evolution 
governed by Eq. \V {C,t)\ attains the critical value, \V\ — 1, i.e., according to Eq. ([6]), the polarization attains its 
critical value, |P| = 1, the further evolution of the system leads to switching from the stable branch of relation ([3]) to 
the unstable one, i.e., from w = — -^/l — |Pp to w = — jf'p- If this happens, the subsequent use of Eq. (|8]) for 

the slowly varying field becomes irrelevant, because the decay of the unstable state may be fast, making it necessary 
to get back to the use of the full system ([T]), ([2]) (which is beyond the scope of this work). 

Asymptotic equations for a single slowly varying amplitude can also be derived near other edges of the two gaps 
that may be occupied by solitons in the full system of Eqs. ([T]) and However, in other cases the eventual equation 
reduces to the usual cubic nonlinear Schrodinger equation (CNLSE). On the other hand, it is relevant to compare 
the NPSE in the form of Eq. ([5]) with the equation which was first derived, under the same name (NPSE), as the 
ID asymptotic form of the Gross-Pitaevskii equation (GPE) for the wave function of a self-attractive Bose-Einstein 
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condensate (BEC) trapped in a cigar-shaped potential In the notation similar to that adopted here, the EEC 
equation (alias the Salasnich equation) is 

Both equations, ([5]) and ((T^ . give rise to a singularity when the local amplitude attains the critical value, \V\ = 1. In 
the case of Eq. (|12p , this singularity leads to the collapse of the wave function, which is a property inherited from the 
full GPE for the self-attractive BEC in the three-dimensional space • The purport of the singularity in Eq. ([8|) is 
demonstrated below: hitting the critical amplitude, the system switches into the unstable phase, which is represented, 
in terms of the underlying RABR model, by square root ([3]) with the upper sign. 

Another type of the NPSE was derived in Ref. [l^ as the effective ID reduction of the GPE for the self-repulsive 
BEC. It seems as Eq. (jS]) with e = -1-1 and -^1 — replaced by y/1 + Of course, such an equation does not 

give rise to bright solitons. Nevertheless, it can generate bright gap solitons, if this nonlinearity is combined with a 
periodic linear potential (the optical lattice) [T9| . 

The rest of the paper is organized as follows. In Section II, we report analytical results for bright and dark solitons 
in Eq. ([S]) with e = +1 and —1, respectively, as well as for cnoidal waves in the former case, and for the switch of the 
system into the unstable phase, in both cases. In the same section we also report results of simulations confirming the 
stability and instability of the bright solitons, as predicted in the analytical form by means of the Vakhitov-Kolokolov 
(VK) criterion. Simulations of 2-soliton states and collisions between moving bright solitons are reported in Section 
III. The paper is concluded by Section IV. 



II. ANALYTICAL RESULTS 



A. Bright solitons 



Solutions to stationary equation (jlip can be represented by means of the formal energy integral, 
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fdW 



const. 



(13) 



For bright solitons, with (|CI = co) ~ 0, which correspond to e = -f 1, as said above, and, accordingly, const = 2/3 
in Eq. (|13p . W{() attains its maximum value (A) at the center of the soliton, where dW/dC, vanishes. Therefore, 
amplitude A can be found by setting dW/dC, = in Eq. (|T3l) with e = +1: 



A' = 




(1-x) 



(14) 



While the bandgap where bright solitons may reside is, formally, semi-infinite in the framework of Eq. ([8]): x < Ij 
the solitons actually exist in a finite interval, 2/3 < x < li in which squared amplitude ([M]) varies from 1 to 0. At 
the limit point of x = 2/3, the soliton solution can be found in an explicit form: 



X=2/3 



(C) 



ICI 



2tanh ( + In (^\/2 + ij j scch |^-|^ +ln fV2 + 1 



The expansion of solution (jl5|) around zero is 



M^x=2/3(C) = 1 - (1/3) C' + (1/6) v/273 Id' + O (C^) 



(15) 



(16) 



As seen from here, a peculiarity of this solution is that, while both W{C) and its first two derivatives are continuous 
at = 0, the third derivative, d^W/dC,^, suffers a discontinuity, jumping from — \/2/3 to -1-^/2/3 as C, crosses zero. 
In this sense, this exact solution may be called a quasi-peakon, usual peakons being solitons with a jump of the first 
derivative at the center [20j . 

In the limit of 1 — x — ^ 0, the smallness of amplitude suggests expanding the radical in Eq. which reduces 
the equation to the usual CNLSE, and the soliton solutions, accordingly, take the following form: 



F = 2yr 



^'^'seeh 



(17) 
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Solitons close to those given by Eqs. ([15]) and (|17|) are displayed in Fig. [T] 

The stability of the solitons can be predicted by means of the Vakhitov-Kolokolov (VK) criterion, dN/dx < 0, 
where N is the norm defined as per Eq. ([9]) [l^l- The numerically calculated curve N{x) for the entire soliton family 
is displayed in Fig. [51 along with the amplitude, A{x), as given by analytical expression ([T3|) (which completely 
coincides with its numerical counterpart), and the soliton's width, L{x), defined as = N^^ ('^W^{()d(. It is 
seen in Fig. ^a) that the VK criterion predicts the stability of the solitons in interval 

Xcr ~ 0.7120 < X < 1, (18) 

and instability in the remaining part of the existence region, 2/3 < x < Xcr- 

Direct simulations of the evolution of perturbed bright solitons, performed in the framework of Eq. ([5]), corroborate 
this prediction: strong perturbations added to VK-stable solitons gradually fade, leaving the soliton intact, as shown 
in Fig. ^a.). On the other hand, the evolution of perturbed VK- unstable solitons does not lead to their immediate 
destruction. Rather, as shown in Fig. [H^b), the amplitude of the soliton grows, crossing the critical level, \V\ = 1. 
Then, formal continuation of the simulations shows a blowup [see, e.g., the panel corresponding to < = 4 in Fig. [31[b)], 
which is a manifestation of the fact that the model as a whole becomes unstable after hitting the critical level, see 
below. 

B. Cnoidal waves 

In addition to the bright solitons, a subfamily of exact periodic solutions in the form of cnoidal waves (expressed 
in terms of the Jacobi's elliptic functions) can also be found from energy equation with e = +1, by setting const 
= X on its right-hand side. First, in the case of 

< X < 2/3, (19) 

the period of the cnoidal solution is defined as an interval of coordinate ^ in which a continuously varying elliptic 
function, sn(^/\/3, fc) with modulus 

k = V(l/2) [1 + (3/2)x], (20) 

takes values that are not too small: 

(V2fc) = (1 + 3x/2)"'/' < sn (C/V3, fc) < 1. (21) 

DJS 
0.7 

6A 
02 
03 
0.1 

%0 -SO -6Q ^0 -^P 3(t 40 SD ?fl 1(11) 

FIG. 1: The solitonic profiles corresponding to x = 0.67 and x = 0.99 (dashed and solid curves, respectively). The numerically 
found solutions are indistinguishable from, severally, the exact "quasi-peakon" sohition (|15p found for x = 2/3, and the 
approximate solution (|17p which is relevant for 1 — x ^ 1- 
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In this interval, the solution is 

W{0 ^ 2k sn (C/V3, fc) dn (c/V3, fc) , (22) 

where dn is the other standard elliptic function, and the entire solution is built as a chain of the so defined periods. 
Note that k given by Eq. ([20| takes values l/-\/2 < fc < 1 if x belongs to region (fTOl). As follows from Eq. ([22|) 
and (j21[) . at junction points between adjacent periods, where the left inequality in Eq. (|21l) turns into the equality, 
the solution attains the critical value, W = 1, and its derivative vanishes, dW/d(^ = 0, hence the matching at the 
junctions is continuous for both W{C) and dW/dC, as well as for d'^W/dQ'^ [according to Eq. (fTTj) . d^W/dC,^ = — x at 
W = 1]. Only the third derivative is discontinuous at the junction, jumping between —^/x s-nd In this sense, 

these cnoidal solutions are built as periodic chains of quasi-peakons, cf. the exact soliton solution ([20|) of the same 
type. Note also that the bright solitons do not exist in the entire interval p^ . i.e., the existence regions of this type 
of the cnoidal waves and solitons are separated. 

Exact solutions for cnoidal waves take a different form at x > 2/3 (recall the solitons exist in the region of 
2/3 < X < 1, i-e., the cnoidal waves may coexist with the solitons in this case, although the cnoidal solutions exists 
also at X > li where the solitons cannot be found). In this case, the solution is built of elliptic functions sn and cn 
with modulus 

k = ^2 [l + (3/2)x]-\ (23) 

cf. Eq. (j20p . Note that expression (j23p is the inverse of (PH)) . and it takes values fc<lforx>2/3. The period of 
the cnoidal solution is now defined by inequalities 

1/V2 < sn (v/(l/6) [1 + (3/2)x]C, k) < 1, (24) 

cf. Eq. (PT|) . and the solution is 

W{C) = 2 sn (v/(l/6) [1 + (3/2)x]C, k) cn (v/(l/6) [1 + (3/2)x]C, k) . (25) 

At junction points between adjacent periods, W again attains the critical value, W = 1, the first derivative vanishes, 
dW/d( = 0, the second derivative is continuous, taking value d^W/d(^ = — x, while the third derivative jumps between 
values ±Y^[cf. the jump of the third derivative in expansion (|16p for the "quasi-peakon" ] . 

Actually, the quasi-peakon solution (|15p corresponds to the limit form of both cnoidal families at the border between 
them, X = 2/3. We stress that the above exact cnoidal-wave solutions, which depend on the single parameter, Xj 
represent only particular cases of a general family of periodic solutions, which depend on two parameters, x and const 
in Eq. (and cannot be expressed in terms of the Jacobi's elliptic functions). 
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C. Dark solitons 



An obvious condition necessary for the stability of dark solitons is the modulational stability of the CW (continuous- 
wave) states, 



CW 



"iXt 



(26) 



with X taking values — 1 < x < 0. A straightforward analysis demonstrates that all CWs (|26|) are indeed stable, if 
e = -1 in Eq. ©. 



Solutions for dark solitons, V (C, t) 



W{C,), which approach CW with frequency x at ^ — >• ±oo, can be 



found as solutions to Eq. (|T3)) with e ~ ~1 and const = x~ (l/3)x^- ^ straightforward analysis demonstrates that, 
in the limit case of x = 0, the solution takes a peculiar form of a compacton- shaped dark soliton: 



^(^^ ^ . sn (C/V3, 1/V2) y2^n2(c7V5^l7V5), at |C| < V^K (l/^/2) 
I 1, at Id > VSK (l/^/2) , 



(27) 



where the modulus of sn is 1/ -\/2, and K (^1/ v^) is the corresponding value of the complete elliptic integral of the first 
kind. Previously, dark compactons were reported in several discrete models [2l[ , but we are not aware of solutions 
similar to that given by Eq. ()27|) in continual models. 



[=209 




FIG. 3: (a) Self-cleaning of a stable soliton with x ~ 0.9, to which a strong random perturbation was added at t = 0. (b) 
The evolution of a slightly perturbed soliton with x = 0.69, which belongs to the VK-unstable subfamily. The crossing of the 
critical level, |V| = 1, implies the transition into the unstable phase. The subsequent blowup indicates the loss of the system's 
stability. 
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D. The switch into the unstable phase 

The crossing of the critical amphtude level, \V\ 1, by an evolving solution can be studied in an analytical form 
too. To this end, a near-critical solution is looked for as 

ViC,t)^l^[v,{C,t)+iv2iC,t)], (28) 

where functions vi and V2 are real and imaginary parts of small perturbations around V ^ 1, i.e., f i 2 ^ l.The 
substitution of this expression into Eq. ^ and a straightforward asymptotic expansion yields the following equations: 

iv2)t~M^^^eV2^, = 0, (29) 
(«i)t + («2)cc = 0. (30) 

Solutions to Eq. ([29| make sense if vi {(, t) does not become negative. 



1. The case of e = —1 

First, we consider this issue for Eq. (|5]) with e ~ —1, which admits stable CW states. In this case, a family of exact 
solutions to asymptotic equations ([29|) and (pO)) can be looked for in the following form: 

(C, t) - (1/2) ~at)\ V2 (C, t) = (1/2) [bt' + c±e + (7/2) C'] , (31) 

with constants [i^h,c^^ . Here i = is defined as the moment of time at which vi vanishes for the first time, i.e., 
the critical value \V\ = 1 is attained, in the framework of the asymptotic approximation. Note that the choice of the 
ansatz for vi in the form of the full square in expression pip guarantees that vi {t, ^) remains positive, as it must be. 
The substitution of the ansatz into Eqs. and ([50)1 gives rise to the following relations: 

b= (1-2/3) a, c = -a^, 7 = {l/3)af3, a = ^2/3(1 - 6/3), (32) 

where /? remains an arbitrary parameter, taking values < /? < 1/6. 

Ansatz (pij) . subject to conditions (|32p . yields an exact solution to Eqs. ([^^1) and (PH)) with e = — 1, provided that 
\/2vi (C, t) is realized, when the ansatz is substituted into Eq. (|29|) . as I3C,'^ — at, but not as \at — (3('^\ (the latter 
expression cannot provide for a solution). An explicit form of -y/l — \V\'^, as given by Eqs. p8)) . ((3T|) . and ([32]) . with 
regard to the above-mentioned realization of y^2vi (C, t), is 

Vl - « ^2vi (C, t) = fie - V2/? (1 - 6/3)t. (33) 

As follows from Eq. ((33)) . a/I — \V\'^ does not vanish (i.e., \V\ < 1 holds) at t < 0, which is the pre-critical stage 
of the evolution. At the critical moment of time, t — 0, ^\ — |yp vanishes at point C = 0. Then, as seen from Eq. 

at t > (at the post- critical stage) -^/l — |yp vanishes at two points, Co(^) ~ ±-\/2 {fi^^ — 6) t. In the instability 
domain between these points, 

- ^2 (/3-^-6)t < C < +^2 (/3-1 -6)^, (34) 

which emerges at i = and expands with time as \/t, expression ()33p . i.e., eventually, square root ([3|) [see Eq. (jG])], 
switches from the stable (lower) branch to the unstable (upper) one. Further evolution of the system is expected to 
be strongly affected by the presence of the instability domain, and should be studied by means of direct simulations 
of the full system of Eqs. ([T]) and ([2|), which is is beyond the scope of this work. 



2. The case of e = +1 

In the case when the bright solitons exist, i.e., e = +1, the critical level, \V\ = 1. is attained at the center of 
quasi-peakon (fT5|) . In this case, it is not possible to find an exact solution to Eqs. ([29|) and pO|) describing the 
crossing of |y| = 1 by a perturbed peakon, unlike the above solution given by Eqs. ([5T|) and p2p . However, taking 
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into regard expansion (|T6l) of the quasi-peakon around its center, and the fact that its frequency is x = 2/3, an 
approximate nonstationary solution can be sought for, at smaU |C| and \t\, as 



(1/V3) Id -at] +0 , 
V2 - (2/3)t+v^t|C| - (1/V2) at' + 0{\C\^te), 



(35) 



where a > is an arbitrary constant, and the square root is reahzed as follows: 

V^-\V? « V2«i iCt) « (v^) ICI - V2at, 



(36) 

cf. Eq. (|33p . As seen from Eq. psp . at t — > —0 this approximate solution describes a small perturbation in the form 
of a cuspon introduced at the central point, C = 0. At t > 0, Eq. (p6| demonstrates that the square root switches 
into the unstable branch in the respective instability domain, |C| < \/3at, which expands linearly with t, cf. Eq. (|34l) . 
As said above, Eq. ^ cannot be used after the switching into the unstable phase. 



III. DYNAMICS OF BRIGHT SOLITONS 



It was mentioned above that, in the limit of 1— x — ?• 0, Eq. ([5]) goes over into the CNLSE. A well-known peculiarity of 
the latter (integrable) equation is the existence of higher-order soliton solutions (breathers), the simplest one, 2-soliton, 



being obtained by the multiplication of a fundamental soliton by 2 at the initial moment |23l |. This circumstance 
suggests to simulate the evolution of so produced double pulses in the framework of Eq. (|8]), if the initial fundamental 
soliton is taken with 1 — x small enough. A typical example is displayed in Fig. Ul which demonstrates that, in this 
case, Eq. ([5]) indeed supports robust breathers of the 2-soliton type, although some emission of radiation is observed 
too. 

Collisions between moving solitons is another natural dynamical problem |2J] . The corresponding simulations were 
carried out, as usual, by taking a pair of far separated identical stable solitons [those with x > Xcr, see Eq. (|18|) ] 
and setting them in motion with velocities ±c by the application of the "kick", i.e., multiplying each soliton by 
exp (±icC/2). Figure [5] demonstrates that, as might be expected, the collisions are quasi-elastic if c exceeds a certain 
minimum velocity, Cmin, while, at c < Cmin, the colliding solitons merge into a single pulse which features irregular 
intrinsic oscillations (as seen from the figure, the quasi-elastic collisions may result in a change of the velocities). 

In fact, the latter outcome is observed only for x > 0.96, i.e., close enough to the CNLSE limit. At x < 0.96, 
the merger produces pulses whose amplitude exceeds the critical level, \V\ = 1 (not shown here in detail), which is 
followed by the switch of the pulse into the unstable phase, as outlined above . 

The velocity separating quasi-elastic and inelastic collisions, Cmin, is shown as a function of x hi Fig. [51 Naturally, 
Cmin increases with the decrease of x, as this implies the growth of the amplitudes of the colliding solitons, and thus 
moving farther from the integrable CNLSE limit, where collisions between solitons are completely clastic. 
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We also simulated collisions between solitons with a phase shift of tt between them. In that case, as might be 
expected, the colliding solitons always demonstrate a quasi-elastic rebound, irrespective of the values of x ^-nd c (not 
shown here). 

Finally, three-soliton collisions (between two moving solitons and a central quiescent one) were considered too. The 
results for them (not shown here) are similar to those presented above for collisions between two solitons: quasi-elastic 
passage of fast solitons, and the merger of slowly moving ones. 



IV. CONCLUSION 



This work aimed to present a new physically relevant variant of the NPSE (nonpolynomial Schrodingcr equation), 
which naturally appears in the model of the RABR (resonantly absorbing Bragg reflector), as the equation for gap 
solitons residing near the edge of the bandgap. The equation features the nonlinearity in the form of the radical term. 
Previously, different forms of the NPSE were derived as models of the EEC dynamics in Refs. and [l^. A full 
family of bright-soliton solutions to the present NPSE was obtained in an implicit form, the ultimate solution being 
the explicitly found "quasi-peakon" . Cnoidal waves in the form of a chain of quasi-peakons were found too. The 
stability of the bright solitons is correctly predicted by the VK criterion, which separates the family into stable and 
unstable parts. The ultimate form of the dark-soliton solution, shaped as a "dark compacton" , was also obtained. 

The instability has a peculiar form in the present model: it occurs when the amplitude of the wave field attains 
the critical value (\V\ =1), which is followed by the switch of the system into an unstable phase (which corresponds 




FIG. 5: (Color online) Outcomes of collisions between identical stable solitons with x ~ 0.96, observed with a gradual increase 
of the collision velocity, c. A transition from the merger to quasi-elastic collisions is observed at c = Cmin ~ 0.07. 
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to the inverted atomic population unsupported by the external field), in terms of the underlying RABR model. The 
passage of the system through the instability threshold was investigated analytically. 

The dynamics of the stable bright solitons was further explored by means of direct simulations. In particular, 
2-soliton breathers were found, and the border between the merger and mutual passage of colliding solitons was 
identified. The merger may lead either to the formation of robust irregularly oscillating pulses, or to the switch into 
the unstable phase. 
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FIG. 6: The velocity separating quasi-elastic collisions and the merger of identical in-phase solitons, versus the intrinsic 
frequency of the soliton, x- The inset zooms the region of small 1 — x- The solid and dashed portions of the curve correspond, 
respectively, to the merger of the colliding solitons into breathers remaining in the stability domain (at \V\ < 1), and to the 
case when the pulses produced by the merger (at c < Cmin) attain values \V\ — 1 and thus switch into the unstable phase. 
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